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We report a polarised neutron study of the magnetic structures and phase trasitions in Gd2Cu04 
in low magnetic fields. These experiments have been complemented by integrated intensity measure- 
ments with unpolarised neutrons in zero field. Polarised neutron flipping ratio measurements have 
been made with magnetic fields H = 0.05, 0.10 and 0.5 T in the temperature range 4-20 K. These 
have enabled us to deduce that the anomalous temperature behaviour of the coherent magnetic scat- 
tering from the Cu sublattice, which shows sharp intensity minima at Tci ^ 18K and Tc2 ~ 8 K, 
is due to cross-overs in the sign of the interaction between strongly coupled, weakly ferromagnetic, 
Cu02 layers. At Tci the coupling changes from ferromagnetic to anti-ferromagnetic and long-range 
order between layers is temporarily lost. Tc2 is the temperature at which the Gd moments order 
and a further reorganisation of the interlayer order takes place. The weak ferromagnetism of the 
CuO layers is found to be due to a small rotation of the Cu moments in the same direction as that 
in which their coordinating oxygen squares rotate in the tetragonal to orthorhombic distortion of 
the crystal structure. Further analysis of the flipping ratio measurements has enabled us to model 
the magnetic structures of the zero-field and the field-induced phases of Gd2Cu04. 

PACS numbers: 75.50.Ee 



I. INTRODUCTION 

The magnetic properties of antiferromagnets that con- 
tain two magnetic species are often very complex due to 
the presence of several competing exchange interactions. 
The series R2CUO4 (R = rare earth element), the parent 
compounds of the electron-doped superconductors, are 
examples of such interesting magnetic systems. Stud- 
ies of Nd2Cu04 and Pr2Cu04i"— have shown that the 
Cu^"*" moments order close to room temperature whereas 
magnetic order in the R'^^ sublattice is only estab- 
lished at very low temperatures. However several spin- 
reorientation transitions which have been observed at in- 
termediate temperatures, notably in Nd2Cu04, demon- 
strate competition between Cu-Cu, R-Cu and R-R ex- 
change. Gd2Cu04 is unique amongst these compounds 
in that although it is as easily doped as other mem- 
bers of the series, it does not, as they do, become 
superconducting.^ Additionally, it shows weak ferromag- 
netism below the Cu ordering temperature T/v (Cu) 
w 285 K— . The large absorption cross-section for ther- 
mal neutrons of natural Gd has limited neutron diffrac- 
tion studies of Gd2Cu04. However such studies have 
been made possible by the availability of ^^^Gd enriched 
Gd2Cu04 single crystals and the magnetic structures of 
the Cu and Gd sub-lattices have been determine d 



A. Crystal and Magnetic Structure 

The crystal structure of Gd2Cu04 at ambient tempera- 
ture is derived from the tetragonal Nd2Cu04 structure 
by a small orthorhombic distortion. ^'^ At Tn — 285 K 
the Cu^"*" moments order in the Nd2Cu04 type of anti- 
ferromagnetic structure with magnetic moments parallel 
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FIG. 1: The layer at z=0 of the orthorhombic structure of 
Gd2Cu04 showing the rotation and distortion of the oxygen 
squares which coordinate the Cu^^ ions. 



to the propagation vector r = (i, i, 0)tetr^- The Gd mo- 
ments order at a much lower temperature T/v = 6.4 K— 
to a structure with zero propagation vector in which op- 
positely oriented moments lie in the a-b plane and are 
related by the centre of symmetry. 

The weak ferromagnetism observed below Tn (Cu) in 
Gd2Cu04 is not compatible with a magnetic structure 
with non-zero propagation vector (vis ^^Otot) but is al- 
lowed in the orthorhombic cell for which the propagation 
vector |^0tct=100oitii is a lattice vector. The conven- 
tional space group of the orthorhombic structure is Cmca, 
however in order to retain the c-axis of the tetragonal 
cell it is convenient to describe it using the space groups 
Acam and Bbcm for the two twins that can be formed 
from the tetragonal phase, aorth = atot ± btct; borth = 
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TABLE I: Conditions on the reflection indices h, k and 
I for nuclear scattering by the orthorhombic structure of 
Gd2Cu04. 



Indices 


h + l 
k + l 


Type 


Description 


hkl 


even 
even 


F 


Fundamental: given by 
the tetragonal structure 


hkl 


even 
odd 


B 


Given by B-face centred 
twin only 


hkl 


odd 
even 


A 


Given by A-face centred 
twin only 


hOl 


even 
odd 


SB 


Space group absence 


Ml 


odd 
even 


SA 


Space group absence 



atct T btot Corth = Ctet- The main effect of the struc- 
tural distortion is to rotate the squares of ions so 
that, as shown in figure [1] these ions no longer lie on the 
[110] and [110] axes of the orthorhombic cell. In the A 
face centred cell (A twin) the Cu layer at z = i is dis- 
placed by with respect to that at 2; = and in the B 
face centred cell (B twin) by ^a. Magnetic domains with 
propagation vector (5^0)tot belong to the A twin and 
those with propagation vector (^-^0)tct to the B twin. 
The squares coordinating Cu^^ ions with antipar- 
allel moments are rotated in opposite directions allowing 
weak ferromagnetism. 

The orthorhombic cell will be used for axes and reflection 
indices throughout the rest of this paper. 
The Bragg reflections from the orthorhombic structure of 
Gd2Cu04 are classifled in table U according to the nuclear 
contribution to their structure factors. 



B. Magnetic phase diagram 



A tentative magnetic phase diagram for Gd2Cu04 based 
on magnetisation and neutron diffraction results^ is 
shown in figure [H The weak ferromagnetic phase (WF) 
is stable from Tn (Cu) down to w 20 K in zero field, 
at which temperature the weak ferromagnetic moment 
drops to zero (AFMl phase). The WF phase can be re- 
tained down to 7 K by application of rather small mag- 
netic fields « 0.2 T at 10 K. Below 6.4 K both the Cu 
and The Gd moments are ordered giving the AFM2 phase 
which undergoes a spin flop transition to the SF phase 
on raising the fleld. The variation with temperature 
of the integrated intensities of the 101 and 102 reflec- 
tions to which the ordered copper moments contribute 
strongly show two anomalous minima which coincide ap- 
proximately with the WF to AFMl and AFMl to AFM2 
transitions.^** The object of the present experiment was 
to study the low field behaviour of the magnetic struc- 
ture using polarised neutron diffraction in order to clarify 
the complex interplay between the order parameters of 
the rare-earth and the transition metal sub-lattices which 
leads to the appearance and subsequent disappearance of 
weak ferromagnetism and to the apparent breakdown of 
order in the copper sub-lattice as the Gd sub-lattice or- 
ders. 



II. EXPERIMENTAL 

The polarised neutron measurements were made using 
the diffractometer D3 which uses a spin polarised neu- 
tron beam from the hot-source of the high fiux reac- 
tor of the ILL Grenoble. The Gd2Cu04 crystal, of size 
w 5 X 8 X 2 mm'^ was a piece of the one used for the mag- 
netic structure determinatio n^ It was mounted in a 
thin tailed cryostat which can pass through the hollowed 
out soft iron poles of the D3 electromagnet. An [010] 
axis was aligned parallel to the field direction. Initially 
the crystal was cooled in zero field and the integrated 
intensities of the 101 and 102 reflections measured as a 
function of temperature between 23 and 2 K. A mag- 
netic field of 0.05 T was then applied which is sufficient 
to ensure 82% neutron polarisation. The flipping ratios 
R of two F type reflections {111} and {113}, one SA type 
{101}, one SB type {102} and two each of the A and B 
types {210}, {212} and {211}, {213} were measured in 
the same temperature range. Most of the measurements 
were repeated in fields of 0.1 and 0.5 T. The integrated 
intensities of the SA and SB reflections 101 and 102 were 
also measured before and after the polarised neutron ex- 
periment using the D9 4-circle diffractometer which is 
also on the ILL hot-source. 



III. RESULTS 



FIG. 2: H vs T Phase diagram of Gd2Cu04i^ 



The temperature variation of the integrated intensities 
measured for the 101 and 102 reflections in the three 
different experiments is illustrated in figure |3] In all 
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three the intensity was found to fall abruptly in a small 
temperature range around Tci and Tc2- Although the 
breaks in the curves occur at the same two temperatures 
the relative intensities at other temperatures are very 
different suggesting that the magnetic intensity in a 
particular reflection depends on the previous history 
of the sample. It was found that that the anomalous 
breaks in intensity can be suppressed by magnetic fields 
of as little as 1 T. 

A polarisation dependence of the scattered intensity is 
expected due to the weak ferromagnetism associated 
with the structural modulation which accompanies the 
Neel transition. This modulation has the same wave- 
vector as the magnetic structure so that the magnetic 
reflections are not entirely of magnetic origin'^'^. In the 
orthorhombic cell the Cu moments are not constrained 
to lie parallel to the crystallographic axes and the weak 
ferromagnetism arises from a small rotation of the 
moments away from the axes. Significant polarisation 
dependence was found in the intensities of all reflections 
except those of types SA and SB for which the nuclear 
structure factors are zero. The variation of the intensity 
asymmetry P = [R — 1)/{R + 1) with temperature, in 
fields of 0.1 and 0.5 T, is shown in figure 2] for a reflection 
of each of the F, A and B types. At the lowest field 
(0.05 T) the asymmetry of the A and B reflections starts 
to increase rapidly above about 12 K. According to the 
phase diagram of referenceii this is the phase boundary 
between the AFMl and WF phases in 0.05 T, suggesting 
that the increasing asymmetry is due the presence of the 
WF phase. The F type reflections 111 and 113 also have 
significant asymmetry at temperatures above Tc2. The 
rapid loss of asymmetry at lower temperature results 
from the rapid increase in the polarisation independent 
intensity due to antiferromagnetic ordering of the Gd 
sub-lattices. 



IV. ANALYSIS 

In a centrosymmetric structure like Gd2Cu04 the cross- 
section can only be polarisation dependent when nuclear 
and magnetic scattering occur at the same wave-vector 
and from the same coherent volume of crystal. In the 
case of the A reflections this magnetic scattering must 
therefore also come from the A twin, and for the B re- 
flections from the B twin although magnetic scattering 
from both twins will contribute to the total intensity in 
the F reflections. For the Cu sub-lattice the magnetic 
propagation vector on orthorhombic axes is either 100 or 
010. Furthermore within the orthorhombic cell the Cu 
ions at (000) and (||0 )are related by an a glide-plane _L 
b for the A twin and a b glide-plane _L a for the B twin. If 
the magnetic moments on the Cu ions at (000) and (^^0) 
are given by the vectors Si and S2 and those on the Gd 
ions at (0,0, izcd) and (i,i,±ZGd) by Ti, T2, T3 and 



T4 respectively then the magnetic structure factors of 
the different types of nuclear reflection are proportional 
to the different combinations of these vectors given in 
table ini Assuming that the moments are always perpen- 
dicular to [001] an estimate of the absolute magnitude 
of the magnetic contributions M to the nuclear reflec- 
tions can be obtained from the intensity asymmetry P. 
If M-M* < IF,/ then 

P = 2{My{kl + kl) - M,k,ky)/P,F„ 

where and My are the components of the real part 
of M parallel to [100] and [010] respectively, kz,ky and 
kz are the direction cosines of the scattering vector, Fn 
is the nuclear structure factor and Pi the polarising ef- 
ficiency. The quantities required for the calculation are 
given in table Hill for the reflections which were measured. 

The 211 reflection is of type A, and comes from the A 
twin only. For this twin in the WF phase, the moments 
are nearly parallel and antiparallel to the propagation 
vector 100 so Si — S2 is parallel to x. The intensity 
asymmetry for 211 is positive above 15 K with the mag- 
netic field in the [010] direction. Then, since kxky/Fn is 
also positive for 211, Si — S2 must be negative to give 
positive P. This shows that for the structure of figure [1] 
in which the square coordinating the Cu^"*" ion at 
000 is rotated clockwise, the major component of this 
Cu^"*" moment must lie in the negative x direction with 
the small positive y component parallel to the applied 
field (see bottom LH panel of figure [SJ . This shows that 
the weak ferromagnetism is due to a small rotation of 
the Cu^"*" moments in the same direction as that of their 
coordinating oxygen squares. 

The magnetic contributions (P_F„/2) to the structure fac- 
tors of the 111 and 113 reflections at 0.5 T and 16 K cor- 
respond to -0.53 and 0.75 /is/cell respectively whereas 
the magnetisation measured by Seaman et al^ amounts 
to only about 2 x 10~^ /ie/cell. The intensity asymme- 
try in these reflections cannot therefore arise from mo- 
ments aligned parallel to the field direction [010] and 
must be due to uncompensated moments in the per- 
pendicular direction [100]. Using the form factors and 
geometric constants given in table IIIIl the moments re- 
quired are 0.47^B/Gd atom and -0.03 //b/Cu. The mag- 
netic contributions to the 210, 211, 212 and 213 reflec- 
tions are 0.92, 0.82, 0.55 and 0.32 /ie/cell respectively. 
They can only be accounted for consistently if the mo- 
ments Si — S2 of both twins, at this temperature and 
field, lie in the [100] direction and the contribution from 
gadolinium is small. These results lead to the model 
for the magnetic structures of the WF, AFMl and SFl 
phases shown in figure [S] Strong anti- ferromagnetic cou- 
pling between the Cu ions in 001 planes combined with 
anisotropy due to the structural distortion leads to the 
formation of weakly ferromagnetic Cu02 layers which at 
temperatures above « 18 K, couple along the c-axis with 
their ferromagnetic moments parallel to one another, to 
give the weakly ferromagnetic WF phase. On lowering 
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FIG. 3: Temperature variation of the integrated intensities of the 101 and 102 magnetic reflections from Gd2Cu04 from 2 to 
22 K (a) Initial experiment on D9, (b) Measurements made on D3, (c) Second experiment on D9 



0.3 

>, 

I 0.2 
>, 

05 
< 

■to.i 

c 

Q) 











1 


1 1 








4 






— • — •— 


— • 















5 10 
Temperature (K) 



15 



0.2 



0.1 



2 1 







▲ . 


* 0.50 
i 0.05 


T 
T 












i 













0.15 
0.10 
0.05 




5 10 15 

Temperature (K) 





^4 






2 


1 1 










i 





















5 10 

Temperature (K) 



15 



FIG. 4: Intensity asymmetry (i? - 1)/(_R + 1) in the F:lll, B:210 and A:211 type reflections in 0.05 and 0.5 T || [010]. 



the temperature below 20 K, exchange interactions due 
to polarisation of the Gd ions lead to a reversal of the sign 
of the interlayer coupling so that between 18 and 16 K 
the long-range interlayer order breaks down and conse- 
quently the coherent magnetic scattering drops. At lower 
temperature anti-ferromagnetic interlayer coupling domi- 
nates and gives the AFMl structure in which the weakly 
ferromagnetic layers are stacked anti-ferromagnetically, 
so that there is no net ferromagnetic moment. The mag- 
netic propagation vector in the A twin changes from 100 
to 010 which is not a reciprocal lattice vector of the A 
cell; the reverse is true for the B twin. Note that the 
absence of magnetic scattering in the {100} reflections in 
zero field at any temperature^ shows that the change in 
interlayer coupling is accompanied by a spin reorientation 
by 90°, so that the moment direction remains parallel to 
the antiferromagnetic propagation vector. 
This model provides an explanation for the anomalous 
temperature dependence of the SA and SB reflections 
of figure [3] The magnetic contribution to the intensity 



of the 101 and 102 reflections from each twin in each 
phase is indicated in table IIVI If the twin populations 
are pA and pb — ^ — Pa then the ratio i?i between the 
mean intensity of 101 in the range 12-15 K to that in 
the range 20-30 K should be « Pa/pb and for the 102 
reflection i?2 ~ Pb/pa- However the products R1R2 were 
found to be significantly different from unity which is 
expected when the Gd sublattices become polarised. If 
the Gd moments have the magnetic space-group of the 
AFMl phase and Ui and U2 are the ratios of the unitary 
magnetic structure factors for Gd and Cu of the 101 and 
102 reflections respectively, then the ratio Ffj_ of the Gd 
to Cu moments can be calculated using 

R1R2 = (1 + F^Uifil + F^U2f « 1 + 2F^(C/i + U2) 

so long as i^^ << 1, in which case pa/pb ~ \/Ri/R2- 
Table |V] shows the twin ratio ij = (pA — Pb)/ {pa + Pb) 
obtained using this relationship. The twin ratios rj and 
magnetisation ratios R^ found in the three experiments 
are significantly different which suggests that the twin 
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TABLE II: Combinations of the magnetic moment vectors of Cu and Gd ions contributing to magnetic scattering in different 
types of reflections in Gd2Cu04. 



Type 




A face-centred twin 




B face-centred twin 




T 


Cu 


Gd 




Cu 


Gd 




F 


100 


Si + S2 


C(Ti +T2 +T3 
+zS'(Ti - T2 + T3 


+ T4) 
-T4) 


Si + S2 


C(Ti + T2 + T3 
+iS{Ti -T2+ T3 


+ T4) 
-T4) 




010 


Si + 82 


C(Ti + T2 + T3 
+2S{Tj - T2 + T3 


+ T4) 
-T4) 


Si + S2 


C(Ti +T2 +T3 
+iS{Ti -T2+ T3 


+ T4) 
-T4) 


A 


100 


Si -82 


C(Ti +T2 -T3 
+iS{Ti - T2 + T3 


-T4) 
-T4) 










010 










Si -82 


C(Ti + T2 - T3 
+iS(Ti -T2- T3 


-T4) 
+ T4) 


B 


100 


Si -82 


C(Ti +T2 -T3 
+iS{Ti - T2 - T3 


-T4) 
+ T4) 










010 










Si -82 


C(Ti + T2 - T3 
+tS{Ti -T2- T3 


-T4) 
+ T4) 



with C = cos 2ttIz and S = sin 2nlz 



TABLE III: Constants relating the components of the mag- 
netic structure factor to the intensity asymmetry P for several 
reflections. 



hklkl+kl Kky F„(mb) SCfcd SSfcd 

1 1 
1 3 
1 
1 1 
1 2 
1 3 



0.5513 
0.7462 
0.8033 
0.8110 
0.8314 
0.8572 



0.4543 
0.2582 
0.4008 
0.3849 
0.3434 
0.2912 



-4.0719 3.4894 



38.9524 
-3.7710 

4.7566 
-6.0364 

4.7566 



3.1567 
2.9723 
2.9371 
2.8349 
2.6753 



■3.7962 
6.5163 
6.2291 
■3.2796 
■2.5969 
5.6540 



6.0438 
0.7310 
0.0000 
5.2213 
5.3881 
0.6343 



TABLE V: Twin ratios 77 and magnetisation ratios calcu- 
lated from the data of figure [3] 



Data from fig 


. m a (d9i) 


b (d3) 


C (d92) 


Ri 


0.73(2) 


1.06(6) 


1.74(3) 


R2 


1.15(2) 


0.980(11) 


0.88(2) 


^/Ri/R2 


0.795(13) 


1.04(3) 


1.40(2) 


V 


-0.114(8) 


0.020(14) 


0.168(6) 


R1R2 - 1 


-0.17(3) 


0.04(6) 


0.54(4) 


Rfj. 


0.042(7) 


-0.010(14) 


-0.135(10) 



TABLE IV: Contributions to the magnetic reflections 101 and 
102 from the WF and AFMl phases of Gd2Cu04. 



Twin 


Phase T S 101 


102 


A 


WF 100 100 pA(4/cS')"g" 
AFMl 010 010 





B 


WF 010 010 
AFMl 100 100 psWcSfq'^ 






with S = |Si| and = (l + {^f) 



populations can change when the Cu^^ ions order. The 
ratio i?^ can be seen to depend on the twin ratio 77 in- 
dicating that the extent to which the Gd moments are 
polarised in the AFMl phase somehow depends on the 
degree of imbalance in the twin populations. 
Figure |6] shows the temperature dependence of the ratio 
between the peak heights of the 101 and 102 reflections 
measured in 0, 0.05 and 0.1 T. The absence of any field 
dependence of this ratio in the WF phase, between 18 
and 25 K, shows that the direction of the magnetic mo- 
ments are coupled sufficiently strongly to the structural 
distortion that the weak ferromagnetic moment of the 



the A twin is not rotated into the [010] direction by a 
field of 0.5 T. The weak hysteresis observed in the mag- 
netisation measurements'''- is then due just to alignment 
of the 180° domains in the B twin. At temperatures 
between 18 and 7 K the energy gain from reorienting 
the weak ferromagnetic moments parallel to a magnetic 
field in the 001 plane may outweigh the inter-layer anti- 
ferromagnetic coupling energy. In this case a gradual 
transition with increasing field, or temperature, takes 
place. With the field applied along [010] the AFMl struc- 
ture of the A twin reverts to the WF structure, but in the 
B twin the moments in alternate planes rotate by 180° to 
give the SFl structure. Si — S2 for both twins now lies in 
the [100] direction and the ratio of the 101/102 intensity 
therefore increases with increasing field as observed. The 
SFl structure proposed would give magnetic scattering 
in the 010 reflection but this cannot be observed in nor- 
mal beam geometry with [010] parallel to the lu axis. 

Below 7 K as the Gd sub-lattices order; there is a rapid 
increase in intensity in both the fundamental and the 
101 and 102 reflections showing that the Gd and Cu 
sub-lattices have the same propagation vector. Since the 
magnetic symmetry of the AFMl structure of figure [51 
Acam' for the A twin, is incompatible with the order- 
ing found for Gd in which Gd atoms at (0, 0, ±z) are 
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FIG. 5: Magnetic structures proposed for the WF, AFMl and 
SFl phases of Gd2Cu04. 
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FIG. 6: (b) Temperature variation of the ratio between the 
101 and 102 peak heights in 0, 0.05 and 0.1 T. 
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FIG. 7: Schematic representation of the structure proposed 
for the AFM2 phase of Gd2Cu04 in the B twin. The Gdl 
atoms are at height i — zod above the layers and the Gd2 
atoms a similar distance below, zod = 0.3393. 



anti-parallel^^, the magnetic space group of the AFM2 
phase cannot have a centre of symmetry at the origin. 
It can be at most Aba2' in which case the Gd sites at 
(0, 0, ±ZGd) are no longer equivalent and the moment 
Ti| need not be equal to IT2I. The polarisation de- 
pendence of the fundamental reflections shows that in a 
field H||[010], Ti + T2 -I- T3 + T4 || [100]. In the AFM2 
phase the magnetic scattering in the F reflections is much 
greater than that in the A and B types, showing that 
Ti - T2 + T3 - T4 » Ti - T2 - T3 + T4 but the 
behaviour of the 101 and 102 intensity ratio below 7 K 
shows that the latter sum Ti - T2 - T3 + T4 || [010] in 
the B twin. All these constraints can be satisfied by the 
structure given in figure [7] for the B twin although nei- 
ther the difference in magnitude of the moments on Gdl 
and Gd2 nor their inclination to the crystal axes can be 
determined precisely from the present data. 



V. SUMMARY AND CONCLUSIONS 

The various magnetic structures which occur in the H-T 
phase diagram of Gd2Gu04 are all built from GUO2 layers 
in which the Cu^^ moments are anti-ferromagnetically 
coupled. In the ideal NdCu204 structure the Cu^"*" ions 
are coordinated by squares of ions lying in the 

planes, but in Gd2Cu04 there is a structural distortion 
in which the squares around alternate Cu^"*" sites 
rotate in opposite directions. The sense of the polarisa- 
tion dependence of the intensity of Bragg reflections with 
both magnetic and nuclear contributions shows that the 
Cu^'^ moments rotate in the same direction as the oxygen 
squares. In the temperature range 285-20 K the weakly 
ferromagnetic Cu02 layers are stacked with their ferro- 
magnetic moments parallel to one another to give the WF 
phase. The interlayer coupling is weak and at lower tem- 
perature increasing polarisation of the Gd ions which lie 
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between the layers favours antiparallel intcrlayer coupling 
and the AFMl structure which has no net ferromagnetic 
moment is stabilised below 18 K. Even a weak magnetic 
field (0.05 T) can reverse this transition in the favourably 
oriented twin and induce a spin flop transition to the SFl 
phase in the other. 

The unusual temperature dependence of the intensity of 
the magnetic reflections in Gd2Cu04 which show minima 
at Tci !v 18 K and at Tc2 ~ 8 K suggests that three- 



dimensional long-range magnetic order is temporarily 
lost at these two temperatures. Two-dimensional mag- 
netic order would be expected to persist and should give 
rise to corresponding rod-like scattering. However, ex- 
perimental measurement of such diffuse neutron scatter- 
ing is not possible with the present Gd2Cu04 crystal 
which has a high absorption cross section even though 
partially enriched with ^^^Gd. 
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